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ABSTRACT 


One of the most important measures of transmission performance is 
the gear tooth root stress (TRS). High root stress can directly damage 
the gear tooth and indirectly affect the life of the gear. This paper 
considered reducing the stresses in the base of the change in the gear 
profile. Accurate calculation of the maximum tooth root stress (TRS) 
and deflection is fundamental to the prediction and optimization of 
gears. The Finite Element Method (FEM) gives accurate results. But, 
the calculation requires a lot of resources and time. Moreover, the 
results obtained from the Lewis, the ISO, and the AGMA methods 
were useful. This paper approved the results of the new method as 
useful as other methods. Therefore, in the current results, a new 
analytical method based on mechanic theory is created using exact 
profile equations to calculate stress and deflection. Firstly, the load 
sharing ratio is considered using the Taguchi method. Finally, tooth 
root stress and deflection are calculated from this profile. The result 
of tooth root stress obtained from the new method is compared with 
the FEM method. The result of the new method is found the 


consistent with FEM method. 
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INTRODUCTION 

Tooth root stress (TRS) is one of the foremost 
performance that indices in gear design and research, 
holds significant importance. A high tooth root stress 
can lead to direct gear damage or have an indirect 
impact on overall gear performance. Accurately and 
ripidly calculating tooth root stress and deformation 
serves as a fundamental basis for gear structure 
design. Currently, there are two primary calculation 
models in use: the FEM model and the analytical 
model. 


The FEM model is a widely adopted approach for 
obtaining accurate and direct measurements of stress 
in the base of gear tooth and deflection. It calculates 
in quantifying the impact of various factors, such as 
different gear periods, shapes, and slope deviations, 
on TRS [1,2]. Furthermore, the FEM model allows 
for the analysis of TRS in gears that cannot be 
effectively assessed by other models, such as thin- 
rimmed spur gears with inclined webs. This model 
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accounts for factors like centrifugal load [3] and web 
angle [4] effects at high speeds. Additionally, the 
FEM model is versatile in reducing TRS by 
optimizing gear parameters to obtain better 
mechanical properties [5]. It is also a valuable tool for 
evaluating the performance enhancement resulting 
from gear tooth profile modifications [6]. 
Furthermore, it can assess the impact of tooth profile 
variations, machining errors, and assembly errors on 
TRS using a three-dimensional FEM model [7]. High 
contact ratio gears present a specific challenge, as 
their load sharing and deformation between teeth 
significantly affect TRS. The FEM model excels. in 
calculating these parameters, a task often difficult to 
achieve with other models. As a result, the scope of 
research encompasses both uniform [8] and profile 
modifications [9] for high contact ratio gears, as well 
as their optimization to reduce TRS [10-12]. To 
improve computational efficiency, new numerical 
models have been developed [13]. 
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This paper analyzed stress on the root of tooth and 
deflection 


generated in anew method. In this article, the stresses 
on the base of the gear tooth and deflection values 
were compared with FEM result. The alternative 
combined effect of the pressure angles, the face width 
of gear, the addendum coefficients of variation, and 
the corner radii adjusted are determined in the new 
model design using the Taguchi method. 


A. Gear Parameters and Material 

The gears data are selected from the lathe machine 
and the mechanical properties of the gear material are 
shown in Table I and Table II. The gear tooth profile 
curve of the new method has used this parameter and 
material properties as the profile shape is shown in 
figure.1. 


TABLE I. Geometric parameters in gear 


parameters 
Parameters Values 
Number of teeth pinion, Tp 40 
Number of teeth gear, Tz 127 
Module, m (mm) 1.25 
Pressure angle, a 20 
Addendum factor, ha 1.25 
Face width, b (mm) 15 


TABLE II. Mechanical properties of gear 
materials 


Parameters Values 


Density of material, kg/m? 7800 
Yield stress (MPa) 986 
Poisson ratio 0.29 
Modulus of elasticity (GPa) | 190 
Brinell Hardness 311 


B. Model of a gear tooth 

The section curve equation can be divided into five 
parts, including the downward curve, downward 
transition curve, upward involute curve, upward 
transition curve, and addendum curve. 


aaddendum circle curve 


upper involutecurve own involute curve 


upper transition curve down] transition curve 


SS ee 


Fig 1. Basic cutter of gear tooth 


The profile equation of the gear tooth is used in 
equation (1) to construct the new model. The section 
of the gear is symmetric about the xg axis, and the 
center of the gear is at the coordinate point O. 


_ mm 
Yow= |g ee 
X,,=h,,m,- X,M,- p; sina 
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Where, Yew is the half length of the cutter pitch line, 
Mn is the normal module, xg is the addendum 
modification coefficient, Xgc, Yge is the center 
locations of edge radius of cutter, rg is the reference 
circle radius, ha is the addendum coefficient, ra is the 
addendum circle radius, @ is the independent variable 
for the curve, p, is theedge radius of cutter, and yg - 


Y ec =Y aw TX gotanat 


Eq: (1) 


angle of the envelope point. 


C. Load distribution ratio 

The first point for the pinion tooth is contact at ‘g,' 
while the second tooth makes contact at 'c.’ This 
results in a sharing of the load between these two 
points, as illustrated in Figure 2, where the applied 
load shares these two gear teeth. 


one -pair of teeth in 
xontact 


Fig 2. Gears contact position 


In a similar manner, points ‘f and 'b' move 
simultaneously with points 'e’ and 'a.' Subsequently, 
‘a’ becomes disconnected, and the full load initially at 
point 'e' is transferred to point 'd' until it reaches point 
'c,, where the newly connected teeth come into 
contact. 


The gear operates with a standard contact ratio, 
meaning that when a single pair of teeth is engaged, 
that specific pair bears the entire load. However, 
when two pairs of teeth are in contact, the transmitted 
load is distributed between the teeth in the 
connection. This load distribution ratio depends on 
both the contact ratio and the stiffness of the teeth at 
the load point position. The load distribution ratio is 
particularly evident at points 'e,''d,' and 'C,' as shown 
in Figure 2. 
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For the cases examined, the load distribution ratios 
remain quite consistent, hovering around R = 0.35, 
and approaching R = 0.65 at the extremes of the 
interval where two pairs of teeth engage. At the very 
limit of the engagement interval for a single pair of 
teeth, the load-sharing ratio reaches R = 1 [15]. The 
load distribution ratio can be expressed as follows: 


RC §)= 0.3540.30° 2 for ban £§ $ En Fal 


(&)=1 for Gntfa-l s€ < & 41 Eq: (2) 
é- Sinn =e 
&— 1 
D. Taguchi Method 
Table.HI outlines the chosen gear parameter levels, 
while it shows the orthogonal arrangement of Lo. In 
this section, the gear profile is modified with respect 
to four parameters: pressure angle, addendum 
modification, edge radius, and face width of the gear 
tooth. However, the face width remains constant. 
Pressure angles are adjusted, and addendum 
modifications vary between 1.25, 1.30, and 1.38, with 
edge radii of fillet at 0.375 mm, 0.5 mm, and 0.625 
mm. 


TABLE III. Level Selection of the process 
parameters 


R( §)= 0.65-0.30 


for FE 71sts € te 


Profile parameter Eexe 


1 2 3 


Normal pressure angle 


19 20 21 
(degree) 
Addendum modification 1.25 | 1.30] 1.38 
coefficients(mm) 


0.375 | 0.50 | 0.625 


A new method has been developed using the Taguchi 
method to consider four different parameter 
variations. The Taguchi method utilizes these four 
parameters as a simple, effective, and systematic 
approach to design optimization, encompassing 
performance, quality, and cost considerations. An Lo 
orthogonal array with 4 columns and 9 rows is 
employed to facilitate the process. Initially, the 
Taguchi method generates nine distinct simulations, 
each utilizing different load distribution ratios and 
corresponding positions. The selection of the 
minimum contact position is preferred to minimize 
applied load[14-15]. Conversely, the maximum 
contact position for a single tooth pair may lead to 
gear failure, as noise and vibration. 


Edge radius (mm) 


Table IV. Taguchi orthogonal Lo array for 
profile level selection 


Pressure Addendum Edge 
No| angle modification | radius 
(degree) (mm) (mm) 
1 19 1.25 0.375 
2 19 1.30 0.500 
3 19 1.38 0.625 
4 20 1.25 0.625 
a 20 1.30 0.375 
6 20 1.38 0.500 
fi 21 1.25 0.500 
8 21 1.30 0.625 
9 21 1.38 0.375 


E. Analytical method of the TRS and deflection 
The highest TRS occurs in the transition curve region 
during single tooth contact, which is in line with the 
gear's working principle. In single-tooth meshing, the 
maximum bending moment typically rises, rather than 
in double-tooth meshing. Furthermore, the involute 
curve region near the root is smaller bending moment 
but transition curve region is larger bending moment. 
The new method introduced in this paper is based on 
these considerations. Figure 3 shows the parameters 
of the TRS and deflection calculation method. 
Following the involute gear meshing rule, the 
direction of the base circle's outer tangent line at the 
meshing point represents the load orientation (F) on 
the gear. The pressure angle at the meshing point 
(Ggr) is the angle between the tangent line and the Yg 
axis, and it varies at different meshing positions. 
Additionally, dr denotes the value of dgiu at the 
meshing point. Taking into account the gear's stress 
characteristics, the standard TRS composed both 
standard bending moment stress and axial 
compressive stress. 


( 


ve (Xgu,Yeu) _ 


Fig 3. Parameters of analytical new model 


The tangent line direction of the main circle outside 
the mesh point aligns with the load direction (F) on 
the gear. The pressure angle at the mesh point (a) 
equals the angle between the tangent line and the ygi 
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axis, and it varies at different positions. Additionally, 
og represents the value at the mesh point. The new 
model was chosen based on modifications to gear 
profile parameters, including adjustments to the 
pressure angle, addendum factor, and edge radius, 
using the Taguchi method. The pressure angle is 
adjusted from 19° to 21°, the addendum modification 
is varied from 1.0 to 1.1 times the module, and the 
edge radius is changed from 0.3 to 0.5 times the 
module. Considering the load characteristics of the 
gears, the standard tooth root stress comprises both 
standard bending moment and axial stress. 


The stress due to bending moment is expressed by 
Equation: 3. [5] 


x _ 3F{ [Xyu (®, )-X ge (®, )]cosa,, ar (9,)}xsina, 


Eq:(3 
2x bx yi(@,) q-( ) 


Where, X¢ is the addendum modification coefficient, a 
is the normal pressure angle, Xgt, ygt is the center 
locations of edge radius of cutter. 


The standard axial compressive stress can be 
expressed by Equation: (4). 


F sina 


z 2xbxy(P, ) Ba: (4) 


The total standard stresses are obtained by using the 
Equations: (5) and (6), these equations are based on 
the base of gear tooth. Total standard compressive 
stress is achieved as follow: 


G,=O,t0, Eq: (5) 
Total standard tensile stress: 
Oy=O_-g Eq: (6) 


F. Deflection of gear tooth 
The deflection of gear tooth is calculated by the 
following eq: (7). 


Deflection, s-= Ea: (7) 


In this equation, 

A=F,xL? 

B=3xEx] 

Where, F; is the tangential applied load (N), L is the 


tooth height, E is the modulus of elasticity (N/mm7), I 
is the moment of inertial of gear tooth (mm*) 


In this approach the strength is modified by various 
factors that produce limiting values of the bending 
stress. 


S,= (0.533H,+88.3) MPa Eq: (8) 


F. Analytical Results 
The load distribution on the gear teeth is determined 
using Eq. (2) at nine different positions. 


1 


0.97 


° 
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Fig 4. Load sharing ratio for all simulations 


Figure 2 illustrates the location of the load 
distribution ratio from the bottom to the tip of the 
gear tooth. The relationship between the load 
distribution along the gear teeth is depicted in Figure 
4. This paper examines the applied load from the base 
to the tip of the gear tooth, considering changes in the 
gear tooth profile facilitated by the Taguchi method. 
The initial contact points are the same, but the final 
contact point differs. Changes in the addendum factor 
and pressure angle result in variations in contact 
lengths and contact ratios. Due to these factor 
adjustments, the maximum load-sharing ratio position 
is identified at 0.24 mm in Simulation 7, represented 
by the blue color. Conversely, the minimum load- 
sharing ratio position is at 0.07 mm in Simulation 3, 
indicated by the red color. Comparing these two 
simulations, the maximum load-sharing ratio position 
decreases by 0.17 mm. The position of the maximum 
load distribution ratio slightly decreases in Simulation 
3. Therefore, Simulation 3 is chosen among the 
others. Gear failure occurs at this position as the 
entire load is applied here. Hence, this position will 
be reduced to prevent gear surface failure. 


The calculation of the maximum tooth root stress 
(TRS) and deflection is performed using the results 
obtained from the load distribution ratio in all 
simulations. Tables V. display the contact positions 
and the maximum stress due to the applied load along 
the tooth flank in simulations conducted with new 
methods. 
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TABLE.V. Tooth root stress (using new method) 
and contact position. 
pointe _— point c 


point g Tooth 


No topoint topoint topoint root stress 

e c a (MPa) 
1 0.80 0.14 0.80 185.84 
2 0.84 0.10 0.84 185.84 
3 0.88 0.07 0.88 185.84 
4 0.75 0.20 0.75 185.84 
> 0.78 0.16 0.78 185.84 
6 0.82 0.12 0.82 185.84 
7 0.69 0.24 0.69 185.84 
8 0.73 0.21 0.73 185.84 
9 0.76 0.18 0.76 185.84 


Table V displays the contact position and maximum 
tooth root stress using the new method along the gear 
tooth. The maximum tooth root stress (185.84 MPa) 
occurs when all loads are applied between points 'e' 
and 'c.' The minimum contact position is observed in 
Simulation 3, making it the position with the smallest 
gear failure due to the applied load. 
200 
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Fig 5. Maximum Tooth Root Stress (MPa) 


Figure 5 illustrates the distribution of tooth root stress 
at various positions, ranging from the tooth root to the 
tip, using four different methods. The stresses exhibit 
a gradual increase from the contact position of 0.5 
mm to 1.42 mm. The stress sharply decreases from 
1.5 mm to 2.42 mm. The maximum tooth root stress 
is observed when the load is applied at 1.5 mm from 
the base of the gear tooth.Although the value of the 
maximum tooth root stress in the new method is 
nearly the same as that in the ISO method, the shape 
of the tooth root stress in the ISO method differs from 
the other methods. 


G. FEM Results 

The FEM method employs ANSYS software. The 
simulation results indicate that the quality of the mesh 
product is influenced by the element quantity, aspect 
ratio, and skewness. Typically, increasing the number 


of nodes (or elements) yields a more accurate result, 
but only up to a certain limit. Consequently, the 
simulation result is governed by the elements and 
nodes, considering the element quantity, aspect ratio, 
and skewness. The values for elements and nodes in 
the FEM method are 1244019 and 292892. 


20.00 60.00 


Fig 6. Use the meshing of pinion 


Figure 6 displays the mesh of the pinion tooth. The 
element quantity for this gear is 0.85652, falling 
within the range of 0.7 to 0.95, indicating good 
quality. Subsequently, the simulation checks the 
aspect ratio, resulting in a value of 1.3809. 


G: New Design A-18 
Static Structural 

Time: 1.5 

4/2/2023 938 AM 


fA) Fixed Support 
B) Force: 1200.N 


Fig 7. Boundary Condition of Pinion 


An aspect ratio of less than 5 is considered good 
quality for most elements (90% and above). The 
simulation also assesses skewness, with a value of 
0.19262 falling within the range of 0 to 0.25, 
indicating excellent quality. Skewness is a crucial 
criterion for mesh quality in a static structure, 
defining the error of cell shape and its equivalent 
volume. A skewness value closer to one suggests 
decreased accuracy. 


The boundary condition is expressed as a fixed 
position in the center of the gear, with an applied load 
of 1200 N acting on the gear tooth as shown in figure 
Te 


Tooth root stress distribution is depicted in Figure 8, 
with the highest concentration observed at the base of 
the gear tooth. 
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Fig 8. Tooth root stress of simulation 3 
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Fig 9. deflection with ANSYS software. 


The maximum deflection, measuring 0.0097 mm, 
occurs at the tip of the gear tooth and expresses in red 
color line. 


a 
180 


H. Conclusion 

The deflections of the spur gear tooth are the same in 
all methods because the maximum applied load value 
is constant. The results of deflection with analytical 
method is 0.0110 mm and FEM method is 0.0097 
mm. The error of these two method is a little different 
(0.0013mm).The gear used in lathe machine was 
analyzed by using analytical method such as new 
method and numerical method such as FEM method 
used ANSYS software. The detail results of tooth root 
stress values are 185.84 MPa for the new model 
method, and 185.22 MPa for the FEM method. The 
allowable strength is 254.06 MPa. So, the maximum 
tooth root stress of the new method is close to the 
FEM method. According to these two results, the 
difference in tooth root stress is 1%. Therefore, the 
tooth root stress of spur gear is sharply decreased than 
the allowable bending stress of the material. The 
analyzed of new profile spur gear will not be failure 
due to these tooth root stresses. 


I. Discussion 
First, the exact profile equation for the involute spur 
gear is established. Then, the maximum tooth root 


stress and deflection of the new method are calculated 
with different tooth normal pressure angles, 
addendum coefficients, and fillet edges, and are 
compared with the FEM method through numerical. 
Finally, the results of the new analytical methods are 
compared to the FEM method. The new method 
utilizes an analytical approach based on mechanics 
theory and possesses stronger theoretical properties. 
When calculating the maximum tooth root stress 
(TRS), both the axial compressive stress and the 
bending stress are considered in the new method. 
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